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ABSTRACT 



This theory permits the calculation of pressure 
distributions over a thin airfoil in steady, inviscid, 
incompressible flow; or given a desired chordwise pressure 
distribution the camber line of an ideal wing can be 
determined. The method treats the circulation about the 
wing -as a continuous vortex sheet of variable strength 
covering the wing planform and trailing downstream to infinity. 
The results of the present numerical solution for the pressure 
distribution solution are not satisfactory. However the 
solution for the camber line of an ideal wing are in reasonable 
agreement with published two dimensional results. 
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I. INTRODUCTION 



This analysis is a continuation of two dimensional 
thin airfoil theory ["Ref. l] and lifting line theory [Ref. 2 ] 
into a three dimensional lifting surface theory, using the 
concept of a continuous vortex sheet over the wing which 
trails off to infinity in the downstream direction. The 
continuous vortex sheet differentiates this method from 
the vortex lattice kernel methods of Multhopp [Ref. 3] , 
Falkner [Ref. 4 ] , Weissinger [Ref. 5] , and Watkins et. al. 

[Ref. 6] and double- lattice finite element methods of 
Lopez and Shen [Ref. ?] and Giesing et. al. [Ref. 8 ] . 

The present method is restricted to steady, inviscid, 
incompressible flow and thin airfoils which are symmetrical 
in semi span. 

The circulation about the wing is treated as generated 
by a continuous vortex sheet of variable strength covering 
the wing planform and trailing downstream to infinity. The 
vortex sheet strength is set by the Kutta condition at the 
trailing edge and requiring no flow through the wing at 
specified control points. 

The flow equations for the wing are resolved into down- 
wash angles, with assumed series solutions for the circulation 
functions . From the known planform geometry and the required 
solution at the control points, no flow through the wing, the 
constants of the circulation functions series are determined. 
The equations are put into matrix form, the matrices being 
formed by numerical summation over the planform, and solved 
by matrix algebra. Once the constants of the circulation 
functions are known all of the airfoil characteristics can 
be calculated. 

This analysis generates the spanwise and chordwise 
pressure distributions over the wing which are necessary in 
designing the wing structure. The solution is also done in 
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reverse, that is the pressure distribution over the wing is 
specified and the shape of the mean camber line of an ideal 
airfoil is determined. 

This analysis uses the consepts of additional and 
reference lift which is of academic interest in under- 
standing airfoil theory. 
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II . General Problem Development 



A. THE WING VORTEX SYSTEM 

The present analysis uses the theory of thin wings to 
describe the characteristics of lifting surfaces in steady, 
inviscid, incompressible flow. The most common planforms 
are those which involve straight leading and trailing edges, 
as illustrated in Fig. 1, but the present analysis applies 
also to wings having curved leading and trailing edges . We 
consider only wings that are symmetrical with respect to the 
x axis. We also restrict the present analysis to lift 
distributions which are symmetrical with respect to the x 
axis . 

The wing and the trailing vortex system associated with it 
are treated as a continuous vortex sheet of variable strength 
y . The sheet strength 'S is regarded as a vector in the 
x,y plane. The wing vortex sheet is treated as lying in or 
very close to the x,y plane and & is everywhere tangent 
to the vortex sheet. 

The sheet strength <T is related to the velocities Vu 
and Vl just above and just below the vortex sheet at a 
given point x,y and to the remote velocity \/oo , as indicated 



in 


the vector diagram, Fig. 2. 


From the vector diagram Fig. 


2. 


Vu = \7eo- i Y* £ 


(1.1) 




Vl = V« - J Y x I 


(1.2) 


It 


then follows that 






( 1/u Vl ) = Vo o 


(1.3) 




(l/u ~Vu) = 


(1.4) 
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or conversely 



y= t* ( Vo - Vu) 



(1.5) 



In Fig. 3, consider the circulation d T around contour 
AOBB 0 A" A . Line segment AOB lies just above the vortex 
sheet, and segment B / Q / A / lies just below the vortex sheet. 



AOS = A' O' &' = ds (1 . 6) 

Then 

dr* f V- ir - 

J C 

d T = l/u'ds + 0 - +0 = ( Vu - Vu ) • as 

df = (f* k ) • ds (I. 7 ) 




It is assumed, which will subsequently be confirmed, that 
the dP in equation (1.7 ) is an exact differential. There- 
fore, there exists a function 

r-- ru>d) 

such that for any small changes in the coordinates 



(t>t) + (yif) (1,8) 

Translating this into vector terms gives : 

dr*(r|i-T|i)-( 7 dx-rjy) = vr-£ (1.9) 

Comparing equations (1.7) and (1.9) shows 

vr = v < £ ( 1 . 10 ) 

or conversely that 

v = kx vr (1.11) 
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Fig. 1 Wing Coordinates 
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i'ig. 2 Velocity Vectors 




rig. 3 Vortex Sheet Element 
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These results confirm that is an exact differential. 

Translating equation (1.11) into Cartesian components 
gives 





(1 . 12 ) 

(1.13) 



The vector function V defined by equation (1.11) or by 
scalar equations (1.12) and (1.13) is non-divergent . 



7 9 V - 







(1 .14) 



B. LIFT 

Newton’s Law applied to a fluid element of fixed mass 

pAXAyA Z. is 



P 



£ A X A y AZ. 



oy 

£>t 



( 2 . 1 ) 



It is assumed that the fluid is non viscous; that is, shearing 
stresses are absent. The resultant pressure force of the 
fluid element arising from the static pressure variation in 
the fluid, is the negative of the pressure gradient multiplied 
by the increment volume. 



resultant = ( - grad p ) 



The pressure force plus the weight is the external force, 
therefore the equation of equilibrium becomes 

f 3 - S rQ ' cf P - ? (2 * 2) 

which is Euler's equation. 

If the fluid is assumed incompressible and equation (2.2) 
is integrated dropping the unsteady term it becomes 
Bernoulli ' s equation 

\/ z 

?4r ’P- 



constant 
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(2.3) 



The pressure difference across the wing at an arbitrary 
point can be found from Bernoulli's equation. 



/ PL- Pu 


\ ~ 1 


/'API 


) Z _ 


. /- 


Vu 




l la. 


/ ~ 1 




^ |/« J 




l/oo 




_ ( Vu + 


Vu) • 


( Vo - Vu ) 


2 V** 


• ( 


y* 


t) 






2. 




v z 




_ 2X ' 


( V 


x k. ) 


2 i • V 


r 




2 




V*. 




14. 






V. 



(2 - 4) 

Equation (2.4) ultimately fixes the lift distribution over the 
wing surface once the circulation function P has been found. 

In the trailing vortex region behind the wing, there is 
no wing surface, and hence no pressure difference across the 
vortex sheet. According to the Kutta condition, the pressure 
difference drops to zero not only behind the wing but also all 
along the trailing edge. This is expressed as follows: 




(2.5) 



This result shows that the vortex lines behind the wing and 
along the trailing edge are all parallel to the x axis, with 
no spanwise components. Spanwise vorticity "jf can occur 

only along the wing surface. 

/ 

In Fig. 4, next page, consider the incremental lift d L 
exerted by the doubly shaded portion of the wing, that lies 
between y and (y+dy) and between x and x. Thus: 

d L = y ? V*, _f ( ) dx d'd 

Xu ^ 

O 

d L - % 14* T dy (2.6) 



This result is a form of the Kutta-Joukowski Law which states 
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y 




chat the force experienced by a body in a uniform stream is 
equal to the product of the fluid density, stream velocity, 
and circulation and has a direction perpendicular to the 
stream velocity. Here P represents the circulation about a 
contour C in the Plane y = constant. Contour C passes in 
front of the wing and pierces the wing surface at point x,y. 

If we now extend the upper limit of integration to the 
trailing x^, holding y constant, we obtain: 
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X X T 

r ( *> y) n; (y) 

clL' , 

dL - e (/» r; dy (2.7) 

thus fr (y) fixes the spanwise lift distribution over the wing 
In order to reduce equation (2.7) to a convenient dimension 
less form, it is useful to introduce the following auxiliary 
notation . 



s = 


total wing 


2 

area, ft 


(2.8) 


2s = 


b = wing 


span, ft 


(2.9) 


c = 


= mean 


chord, ft 


(2.10) 


Ak = 


II 

”ji|^ 

II 


aspect ratio 


(2.11) 


n = 


M 

— = spanwise coordinate 


(2.12) 


ii 


dy 

s 




(2.13) 




d 1— 


section lift coefficient 


(2.14) 




o^«.c 
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J 

<U 


L 

-S 


wing lift coefficient 


(2.15) 



By utilizing the above nomemclature, it is easy to reduce the 
spanwise lift distribution as given in equation (2.7) to 
either of the following two equivalent forms. 

(2.16) 

To find the total lift we integrate over the span, For 
symmetrical lift distribution it is only necessary to 
integrate over the semi-span and multiply by two. Thus: 
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(2.17) 



o O 

C. INDUCED VELOCITY 

The velocity cW P ' induced at an arbitrary point P, 

coordinates x ,y on the wing surface by an element of the 
P P 

vortex sheet which lies at the point whose coordinates are 
x,y is found as follows. The vortex sheet element has length 
ds, width dn and strength Y . The circulation strength of 
the vortex filament involved is = Ydn. 

The Biot Savart Law for this case may be written in the 
form : 



dvv p 



Y x Y dn ds 
4 tt r 3 



(3.1) 



where 



T = i (x-x ? ) + -j ( 4 - HpW k C z -" z p) 



(3.2) 



dn ds = dS = element of the wing area 
Y= (. Yk + i k Yz. - k * V T 



(3.3) 



Neglecting the terms in equation (3.2) and k. V"z. in 

equation (3.3). 






also 



r 4 



i 



i 



k. 



(x-xp) I'q-'dr) 

- — vr 

& 'j 



o 

o 



- k [(*- x ^(^r) + ( 1 



(3.4) 



(3.5) 
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Fig. 5 Induced velocity in the x-y plane 
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Upon substituting expressions (3.2), (3.3), and (3.5) into 
(3.1), and integrating over the entire wing surface S 
( but not yet over the vortex sheet trailing behind the wing) 
the downwash due to the wing vortex is obtained. 



W P = 



!e 

TF 



' £(*-*») TV ^ 3 

s [(*-*.r* ] iA 






(3.6) 



The incremental downwash velocity dwp induced at point 
P due to the semi-infinite trailing vortex filament of strength 
cl Hr which leaves the wing at the point whose coordinates are 
x T and y as shown in Fig. 5 is by the Biot Savart Laws 





y 



( 

4TT ' 



y -yp 



(X-Xp) 



* k 



(X-X & ) 



-$> 



(_L_\ , 







(3.7) 




the 


wing gives 


• 

1 


H 




(3.8) 



The total induced velocity at point P is the sum 



W P = M/p -+■ Wp (3.9) 

where Wp and vvp are evaluated respectively from equations 
(3.6) and (3.7). The first of these integrals represents the 
velocity induced by the vorticity over- the actual wing surface, 
and the second integral represents the velocity induced by 
the vorticity trailing behind the wing. 

The downwash integrals are now non-dimensionalized and 
the Cartesian coordinates x and y are transformed into 
dimensionless relative coordinates f and £ defined as 
follows s 
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= relative chordwise coordinate 



(3.10) 



X-Xt_ 

C- 

d n . 

j- = relative spanwise coordinate 

The following auxiliary definitions are used. 

c = -j— = average chord 

/q Zs. 

™ = — — = aspect ratio 

(?) = — r— = leading edge function 

2! r (7) = = trailing edge function 

1 c 

= C (7) = relative chord function 

C<.?) = X t (T) -2.J?) 




(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 



Note that if the parameter is specified and if any two 
of the three functions involved in equation (3.17) are 
specified, this information completely defines the wing 
planform. 

The circulation function T is non-dimensionalized as 
follows. The number four is inserted into these definitions, 
somewhat arbitrarily, because this reduces the final downwash 
integrals to a particularly convenient numerical form. 

= dimensionless circulation 
function along the wing 
surface (3.18) 

= dimensionless circulation 
function at and behind the 
trailing edge (3.19) 



r *(}>•?) - 






r.d) = 



Tt 

”4 Vco-S 
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The induced velocities Wp and vv r are non-dimensionalized 
to corresponding downwash angles c<^ and «: P " as shown below. 
Those angles are defined as positive if down from free stream 
velocity . 



/ 

0 ** p 



ft 

(Xp 



Wp « la. 
V'eo 



downwash angle point P induced 
by vorticity distribution over the 
wing (3.20) 



— ft . 

Wo « fe. 
V'oo 



downwash angle at point P induced 

by vorticity distribution behind 

the wing (3.21) 



p = total downwash angle induced at point 
P by the complete vortex sheet system 
of the wing (3.22) 



Some useful geometrical relationships are : 



(X-Xp) = s, JL[(X l -Q) -(^up +Cp?p) ~] (3.23) 

(oj-^p') - s ( ^ - *? p ) (3.24) 

Also at any point on the wing 

r(x^) ^ "M/oos r* ( ? > y) (3.25) 

Differentiating equation (3.25) 






Differentiating equations (3.24) and (3.25) 
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= s, -^r [_C A l + (X!l + C f ) d *X " ° 3 



c{x - - ^ 
dy = s. d >1 



(3.27) 

(3.28) 



The prime marks in equation (3.27) denotes differentiation 
with respect to y. Substituting equations (3.27) and (3.28) 
into equation (3.26) gives 



= (t£) 4 + C 1 






n 



c> >1 



ft 



(3.29) 



The coefficients of dz and dy must be separately equal on 
both sides of equation (3.29) since this relation must be 
satisfied for arbritrary values of dz and dy. Hence 



(£)* c = 


) 




(3.30) 




/ 


) - 


(3.31) 


^ r 

Solving those two equations for 


V P 

and - — gives 
^ & 






> 


M 


(3.32) 


(&)- 




(3.33) 


The element of wing area in relative coordinates 


may be 


expressed in the following manner. 






dS = 


>*- s 






ds-- Zs z ( z \ )(£ ) JLdt = 




(3.34) 
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dS 







Fig. 6 Element of Wing Area 



The various quantities of equations (3.6) and (3.8) are 
transformed into non-dimensional notation. 



f)1 1 



= ^(£[(ZsC{)-(X^<ri P )] {{j 

+ h L- ^ ) -C^vC- )] 



^ 4!4 .s« |-{[(.s:vC£)-(X u .-CpIp)-C«>>)«^!)](|f- 

- W-'M C(x^-) dfc!* (3.36) 



(x T -x») = s. ^ [r- - ^ - c P fh 



(3.37) 
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Study of the algebraic structure of these relations 
shows that the final downwash integrals can be reduced to a 
relatively concise form by adopting the following auxiliary 
nomenclature. Let 



[~z[(X l +Q.[)-(^p<f!p) 1 Z + (l-tlr)* J* z 



F* = 






(£-*■ C f | p)] Z + (r(->tr>) Z J 



- 0 - 












(3.39) 









(3.40) 



The downwash angle integrals equations (3.6) and (3.8) can 
now be expressed in the following form. 




(3.41) 



(3.42) 



D. THE CONDITION OF NO FLOW THROUGH THE WING 

The shape of the mean wing surface near the x,y plane 
may be expressed by the function 



z = z(x,y) 



(4.1) 



A unit vector nV normal to the above surface at the arbitrary 
point P on that surface, in terms of its Cartesian components, 
is given by the expression 
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_ -( jit )? i. - (.4^i i *l!lk 

n " = (4 

Since and (- 5 ^-)* are negligible in comparison with' 

unity, equation (4.2) is adequately approximated by the 
linearized version 



ho - 



= - ^ 



z 



0 ) 



(4.3) 



The resultant velocity vector at point P may be closely 
approximated by the expression 

V P - T V'ok Cos.<* •+• k. sim X - )e tA/p (4.4) 

Equation (4.4) implies that the induced velocity w/p is taken 
as positive if down. 

Since x is a small angle equation (4.4) may be linearized 
to 

Vr - <- V<*> + k (l^cc-wV'} (4.5) 



The condition of no flow through the wing now requires that 






. n - o 



(4.6) 



thus 



( 4 - 7) 

()4,-<-wp)= o (4.8) 



hence 
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r.lri - 

^ 'bx - 



<K - 



Wp 



(4.9) 



but 




p - o( 



p +- ot. 



therefore 



(4.10) 




The overall wing equation is now found by substituting 
the integrals of equations (3.41) and (3.42) into equation 
(4.11). The result is the basic wing equation. 




o< 





(4.12) 



The influence functions F^ , F^, and F^ are as defined by 
equations (3.38), (3.39) and (3.40). These three functions 
depend only on the planform and aspect ratio of the wing. 
Theoretically, equation (4.12) must be satisfied at all points 
in the wing surface for an exact solution. However, an 
adequate approximate solution can frequently be obtained 
by satisfying equation (4.12) at a sufficient number of 
discrete control points suitably distributed over the wing 
surface . 

E. REFERENCE AND ADDITIONAL LIFT DISTRIBUTIONS 

Let the wing slope function be known at every 

point , tj p of the wing. Also let the entire leading 
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edge of the wing lie exactly in the x-y plane and the 
trailing edge lie in the x-y plane at midspan. In general, 
however, z is not necessarily zero at other points along 
the trailing edge. Let 17* (t,?) be a function which 
satisfies equation (4.12) when has a particular value . 
The subscript r stands for reference condition. Thus at the 
reference condition equation (4.13) becomes 




The general solution or e Q ua tion (4.13) can now be 

written as the sum of the reference solution and an additional 
solution in the form 

ni.'?) = (5.2) 



To demonstrate that equation (5.2) is valid, and to obtain 
the equation that governs the form of the additional solution 
function we proceed as follows. Substitute equation (5.2) 
into equation (4.12). Subs tract equation (5.1) from the 
result. Finally divide through by the factor (^oc- cc r ). The 
result is 



i i 



O = 



i 

TT 



r 



Z l ^ ^ \ 

F. \Tf) 



( ^ > 









+ 



r 



j 

-i 






cj fr 



V 

To- 



d v\_ 



) Jf 



(5.3) 



Thus equation (5.1) governs the reference solution or the 
particular solution while equation (5.3) governs the 
additional solution. The specific value of in equation 

(5.1) is fixed by the particular condition that 
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z = 0 at |= + 1,^=0 



(5.4) 



given that 



z = 0 at £ =0 for all f ( 5 . 5 ) 

The reference solution depends not only on the wing 
planform and aspect ratio but also on the wing slope function 
and the angle of attack . However, the additional lift 
function involves only the wing planform and aspect ratio. 

It will be shown later that the wing lift curve slope 
depends only on the additional lift function n*T and is 
independent of the reference solution Tp • Therefore, 
the wing lift curve slope depends on wing planform and 
aspect ratio only, and is independent of the wing slope 
function. 

To specify the design of a wing it is necessary to fix 
two planform functions and c(l) » the aspect ratio 

AK , and the wing slope function (- 57 -) p at all points 

r[ P on the wing. For a wing of specified design, equation 
(5.1) must be solved for the initially unknown function 
fV and equation (5.3) must be solved for the initially 

unknown function pp (1 j?) • 

In the case of equation (5.1) the above procedure can 
also be partially reversed, that is, the function pp {]>*[) 
can be specified arbitrarily, within certain limits, and the 
equation can then be solved for as the unknown. This 

reversed solution procedure had two advantages. Firstly, it 
means that the lift distribution over the wing in the 
reference condition can be stipulated independently and that 
the wing slope function can always be found such as to yield 
the desired reference lift distribution. Secondly, it 
happens that equation (5.1) is much easier to solve when the 
wing slope function is the unknown than when the reference 
lift distribution is the unknown because, in the latter case, 
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the unknown function is under the integral sign. If 



is the unknown, equation (5.1) can be solved by direct 
integration; usually this must be performed numerically. 
Repeated integrations of equation (5.1) are required, a 
separate numerical integration being needed for each discrete 
point at which the wing slope function is being evaluated. 

In order to represent the wing slope function adequately 
it is necessary to choose a sufficient number of discrete 
points suitably distributed over the wing. 

Unfortunately, the above reversed solution procedure 
cannot be applied to the additional lift distribution ToT 
as governed by equation (5.3). Of course, the equation can 
still be solved, at least approximately. 

The choice of the reference angle of attack <x r which 
fixes the corresponding reference solution Vr is to a certain 
extent arbitrary, The most common choice is to set c<v equal 
to the angle of zero lift of the wing, that is to the angle 
which yields a wing lift coefficient of zero. For a 
symmetrical lift distribution the lift coefficient is 



The restriction that the lift coefficient equal zero does 
not necessarily mean that the integrand is identically zero 
all along the span, just that the net area inder the curve 
sums to zero. Thus at = 0, portions of the span may be 
exerting upward forces, while other portions are exerting 
downward forces . A wing that behaves in this way at = 0 
may be said to be aerodynamically twisted. Conversely a 
wing that at = 0 yields c^-^- identically zero for all 
values of ^ may be said to be aerodynamically untwisted; 
this however does not imply the absence of camber. 

The reference condition for this analysis is that which 
corresponds to the angle of attack for which the wing lift 
coefficient equals zero. This particular lift distribution 
is also commonly labeled as the basic lift distribution. 




(5.6) 



O 
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F. ANGULAR COORDINATES 



A 

?>< 

made 



then 



further shift from the linear relative coordinates 
to corresponding angular coordinates 0 , e , is now 

See Fig. 7. The basic conversion relations are 



2 - 


h ( i 


- COS P ) 


(6.1) 


f = 


cos e 




(6.2) 


oil 


h sin 


/ 

p 


(6.3) 


r 

ii _ 

d© 


- sin 


© 


(6.4) 


n 


) = 


2 ^ ( \ 


(6.5) 


— ( 


) = ' 


' — ( 'l 

Si^sl 0 ^ 3 


(6.6) 




Fig. / Angular Coordinates 
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The downwash integrals can now be converted as follows 
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O TT 
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f t s.m r 2 (HI) 
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o o 



TT IT 
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O o 



-«r=w[ F >(0H- 
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where 



_L 

tt 






* (^) ^ 



( 6 . 7 ) 



(• 

( 6 - 8 ) 



Gi = sin 6 



G 0 = - _t£ sin <2^ 
2 X 



G 3 " " F 3 



( 6 . 9 ) 

( 6 . 10 ) 
( 6 . 11 ) 



Now by refering to the above relationship and to the previous 
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equations which define the functions , and F^ it is 

possible to summarize the final wing equations as the 
culmination of a series of systematic calculation steps. 
The following scheme of auxiliary notation is used. 





l = 


h ( 1 - cos <f> ) 




(6.12) 




Ip = 


hi 1 - COS ) 




(6.13) 




7 = 


cos 




(6.14) 




= 


cos ©p 




(6.15) 




(it = 


iu( e V C C^ I 




(6.16) 




hip = 


X u (_ ©p) + C 




(6.17) 




(u = 


■4 L z >- - 1 Lc- 0 




(6.18) 


Depending on planform, the derivatives 


d © 


and i£. 
d © 


may be multi-valued at midspan, however 


for a 


symmetrical 


planform they are set to zero. 






h 3 = 2 


sin 


Hip ) - hz j 




(6.19) 
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(ht- hip ) J -f {/(- ^p) 


[ 


(6.20) 
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(6.21) 


h 5 = . 






(6.22) 
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(6.25) 
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Os - — 



( 6 . 26 ) 




The basic wing equations become the following 



Reference lift 



IT tt 

r 



= ° <r + if*. f |L G ' 1 ^ 



■ Te i j 



© 



4 - 




( 6 . 27 ) 



Additional lift 



O - 








c> 



G. APPROXIMATE SERIES SOLUTION 
The reference and additional 
represented by double series of 



( 6 . 28 ) 

distributions can be adequately 
the following form. 



W £ 



f r ( <^ ) © ) - Sr s 



Oro 

ftl S 






( 7 . 1 ) 



( 7 . 2 ) 



For a spanwise symmetrical lift distribution, only odd values 
of s are used. Ref. 1. There are (R-rl ) distinct functions Hr 
and ( ^ ) distinct sine functions in the above series. 

The total number of initially undetermined constants in each 
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of the above series is therefore 

T = ( R +- 1 )•( ) (7.3) 

If there are (R+l ) chordwise stations and ( ■ S> 'gy - ) s panwise 
stations over the semi-span, this establishes T control 
points. By satisfying each governing equation at each of the 
T control points we can solve for each set of T constants . 
Once the constants Brs and Ars are found, all aspects of the 
wing performance can be readily calculated including lift 
distributions, slope of the lift curve, induced drag, etc. 

The character of the H functions is established as 
follows 



AP _ 2 ( 

Vco K ** ) 








(7.4) 


>> r zLfi^v<*> (ir*\ 

CL v > 








(7.5) 


therefore 










AP t 4-fc\['hV i ’\ ( 4-/^ ' 


\ 2. 


( >r *) 




(7.6) 


' CL ^ C_ 


/ Sirvi 9$ 






For the additional lift, equation 


(7.2) 


is substituted 


into 


equation (7.6) 










1 (Vo S ••• 


r 2 




1 S/M S © 






M0 H 


(7.7) 


For the reference lift, equation 


(7.1) 


is substituted 


into 


equation (7.6) 










R+i -s 

(f-)r=0 H Br ^ 

f-O \ ) " 


z 




SiM s © 




>(Ni <$ 


l\ 


(7.8) 



It is now possible to make use of the results of thin 



airfoil theory and set up the functions in equations ( 
and (7.8) in an analogous form. Ref. 2. Thus for the 
case of r = 0 we use 

4 __ 2 . _ \ 4- cos <& 

'b | d ~ S/M jz5 

and for r ^ 0 we use 




Equations (7.9) and (7.10) can be integrated to yield 
For r =0 



Ho(^ = 4r ( + S/M J* 'l 

For r = 1 

H, 

For r > 2 

|_j r = JL f siNj(r-i) _ sim (r-vi) gS 7 
r ^ l r - 1 r+ 1 i 

At the trailing edge, g6 = tt , these relations yield 

Ho Lrr)~ TL 

HAtt) = ~ 

Hr-(TT) - O 



.7) 

special 

(7.9) 

(7.10) 

(7.11) 

(7.12) 

(7.13) 

(7.14) 

(7.15) 

(7.16) 



By substituting these values into equations (7.1) and (7.2) 
the corresponding expressions for the trailing vorticity 
are obtained. 



(7.18) 



fn-a. L & ) ~ ( 2- Aoi + S "J 5© 

An ideal wing is one for which = 0 at 

the leading edge. This requires that Bos = 0 for all s. 

If the planform function is arbitrarily- prescribed the wing 
should not be presumed to be an ideal wing and all terms in 
Bos should be included in the solution. The additional 
lift is always characterized by the singularity (•^2)— 00 
at the leading edge so that the terms in Aos must always 
be retained. 
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Ill DESCRIPTION OF NUMERICAL PROBLEM SOLUTION 



A. CALCULATION SCHEME 

The calculation scheme in this analysis uses two distinct 
meshes. There is a coarse mesh of control points, and a 
fine mesh for the purposes of accurate integration. It is 
desirable with respect to both meshes that the ratio of the 
number of spanwise stations over the semispan to the number 
of chordwise stations should be roughly equal to the aspect 
ratio divided by two. This gives roughly equal resolution 
in the spanwise and chordwise directions. 

The total angular interval, 0 to n - in each case is 
subdivided into equal sub-intervals and the calculation point 
located at the center of each sub-interval. The calculation 
points associated with the fine mesh must not coincide with 
any of the control points associated with the coarse mesh 
because this would cause the functions and F 2 to assume 
an indeterminate form. While this indeterminacy can be 
resolved it is best avoided. A good way to avoid this 
difficulty is to make the fine mesh an even submultiple of 
that of the coarse mesh. Then the control points will 
always lie along intersections of the boundaries of the area 
elements which comprise the fine mesh. See Fig. 8. 

The double integration over the wing area reduces in this 
scheme to a summation over the points of the fine mesh. The 
number of terms in the series being equal to the number of 
control points. 

If the pressure distribution is to be found for a given 
wing the coefficients of the reference and additional lift 
series are solved for first and then the pressure distribution 
can be calculated. 

If an ideal wing camber is to be defined from a given 
chordwise pressure distribution, the wing equation is solved 
for (|rr)p • The terms in M* are dropped which eliminates 
the leading edge singularity, and an equal number of terms is 
added on to the end of the series. In other words, Bos is 
set equal to zero. The induced drag of the wing can be 
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INTEGRATION POINT 




Fig. 8 Calculation Meshes 
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reduced to an absolute minimum if the terms in sin s e> 
are eliminated retaining only the term in sin . This 
defines an elliptical spanwise load distribution. These are 
desirable features to incorporate into the design of an 
ideal wing, and it also simplifies the numerical solution. 

B. MATRIX FORMAT 

By analogy with equations (7.7) and (7.8) the reference 
and additional lift distributions may be written in matrix 
format . 





■l - l>] H 


(8.1) 




J [a] 


,(8.2) 



The matrix format of the basic wing equations (6.27) and (6.28) 
become 



("If ~* r j - (8 - 3) 

[i]c-[r][a] (8.4) 

The matricies TqI and in the above equations are known 

T by T matricies whose elements depend only on wing planform 
and aspect ratio. 

Solving equations (8.3) and (8.4) for the initially unknown 
constants B and A gives by matrix inversion 

— j 



(8.5) 

( 8 . 6 ) 
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Equation (8.6) shows that the A constants depend only on 
wing planform and aspect ratio while the B constants in 
equation (8.5) also depend on the wing slope function. 

The points used in the Q matrix need not be the same as 
the control points used in the R matrix, but may be any 
point on the wing. In this analysis the points for the Q 
matrix were picked as six stations across the semispan with 
eighteen chordwise stations at each spanwise station. This 
permits a good representation of the functions T and (^“) ^ _ 

If the reference pressure distribution over the wing is 
specified, coefficients of the lift distribution series can 
be found, by equation (8„7)o Then the wing equation (8.8) 
can be solved by direct numerical integration. 






(8.7) 
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IV RESULTS 



A. SOLVING FOR THE WING PRESSURE DISTRIBUTION 

The numerical results of the pressure distribution solution 
were not satisfactory. The pressure distributions and circu- 
lation functions calculated by the computer program were not 
smooth or in accordance with [Ref. l] and [Ref. 10] .It is 
believed this is a result of difficulties with the first 
chordwise term, Bos. The term by term downwash is plotted in 
Fig. 9 , at a spanwise station near midspan for a wing of 
elliptical planform and an aspect ratio of 20 using 5 chord- 
wise stations and 10 stations over the span. The integration 

mesh was four times finer than the control point mesh. The 
first term should uniform downwash. the rest of the terms 

are of the general nature expected. 
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B. SOLVING FOR THE WING CAMBER 

The numerical results of calculating the camber of an ideal 
wing given the desired chordwise pressure distribution were 
satisfactory. Fig. 10 compares the results of — x r for a 

wing of aspect ratio 20 and rectangular planform using a 15 
term lift distribution series at 5 chordwise points near mid- 
span, with the two dimensional results of [Ref. 10.] It is 
felt the differences at the leading edge can be resolved by 
using more terms in the lift distributions series and more 
integration points. The depicted pressure distribution was 
imposed at all spanwise stations. The results are normalized 
to a lift coefficient of one. Table 1 shows the results at 

various spanwise stations. 

Table 1 

AR=20 15 term lift distribution series rectangular planform 
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0.1400 


0.1669 
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-0.0967 
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-0.1 884 


-0.2142 


-0.2256 




Fig. 10 Camber Line 
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V. CONCLUSIONS 



The numerical implementation of the theory needs 
additional work to achieve the correct computer solutions 
for the pressure distribution solution. The areas of 
greatest interest are the assumed series solutions for 
the additional and reference lift, especially the first 
chordwise term, Bos. The camber line solution is consistent 
with two dimensional theory near midspan of a high aspect 
ratio wing and allows the camber line to be determined at 
various spanwise stations for any aspect ratio. It is 
felt that this theory has merit and deserves further study 
to correct the numerical solution of the pressure 
distribution . 
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APPENDIX A. DERIVATION OF SUBROUTINES 



A. SUBROUTINES FOR THE PRESSURE DISTRIBUTION SOLUTION 

QMAT calcualtes the Q matrix of equations (8.1) and (8.2) 



s 






(7.7) 









(A . 1 ) 
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(A. 2) 



Si = S/aJ s <=> 
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c 



(A. 4) 



RMAT calcualtes the R matrix of equations (8.3) and (8.4) 
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C 5 

B. SUBROUTINES FOR THE CAMBER LINE SOLUTION 

QMAT calculates the Q matrix of equation (8.1) which has 
been simplified with a single term in theta, an elliptical 
spanwise lift distribution, and the first chordwise term is 
dropped . 
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(A. 16) 



<p( r) = . H H 1 ♦ Si 

c. 

CAMBER is a numerical integration of equation (6.27) 
simplified for an ideal wing. 
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